Orbital fluctuations in the diff'erent phases of LaVOa and YVO3 
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We investigate the importance of quantum orbital fluctuations in the orthorhombic and monochnic 
phases of the Mott insulators LaVOa and YVO3. First, we construct ab-initio material-specific t2g 
Hubbard models. Then, by using dynamical mean-field theory, we calculate the spectral matrix as a 
function of temperature. Our Hubbard bands and Mott gaps are in very good agreement with spec- 
troscopy. We show that in orthorhombic LaVOa, quantum orbital fluctuations are strong and that 
they are suppressed only in the monoclinic 140 K phase. In YVO3 the suppression happens already 
at 300 K. We show that Jahn- Teller and GdFeOs-type distortions are both crucial in determining 
the type of orbital and magnetic order in the low temperature phases. 
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The Mott insulating perovskites LaVOa and YVO3 
exhibit an unusual series of structural and magnetic 
phase transitions (Fig.[T]) with temperature- induced mag- 
netization reversal phenomena f]\ and other exotic prop- 
erties 0, [3| . While it is now recognized that the V- t2g 
orbital degrees of freedom and the strong Coulomb re- 
pulsion are the key ingredients, it is still controversial 
whether classical (orbital order) [E 0, 0, d, 0, Hi or quan- 
tum (orbital fluctuations) [3, 0] effects are responsible for 
the rich physics of these vanadates. 

At 300 K, LaVOs and YVO3 are orthorhombic param- 
agnetic Mott insulators. Their structure (Fig. [2]) can be 
derived from the cubic perovskite ABO3, with A=La,Y 
and B=V, by tilting the VOg octahedra in alternating 
directions around the b-axis and rotating them around 
the c-axis. This GdPeOa-type distortion is driven by AO 
covalency which pulls a given O atom closer to one of its 
four nearest A-neig hbors Since the Y Ad level is 

closer to the O 2p level than the La 5d level, the AO co- 
valency increases when going from LaVOs to YVO3 and, 
hence, the shortest AO distance decreases from being 14 
to being 20 % shorter than the average, while the angle of 
tilt increases from 12 to 18*^, and that of rotation from 7 
to 13° [mil. Finally, the A-cube is deformed such that 
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one or two of the ABA body-diagonals is smaller than 
the average by, respectively, 4 and 10 % in LaVOa and 
YVO3. These 300 K structures are determined mainly 
by the strong covalent interactions between O 2p and the 
empty B Cg and Ad orbitals, hardly by the weak interac- 
tions involving Bt2g orbitals, and are thus ver y s imilar 
to the structures of the t^g La and Y titanates jiol . ITU . 

The vanadates, however, have a much richer phase 
diagram than the t2g titanates. At, respectively, 140 K 
and 200 K, LaV03 and YVO3 transform to a monoclinic 
structure in which c is turned slightly around a whereby 
the two subcells along c, related by a horizontal mirror 
plane in the orthorhombic structure, become independent 
(Fig.Hl). Most important: a sizable (3-4%) Jahn- Teller 
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FIG. 1: Temperature-dependent structural and magnetic 
phases of LaVOa and YVO3. The lines show LDA-I-DMFT 
(quantum Monte Carlo [2l|]) results for the occupations, n, of 
the three t2g crystal-field orbitals, 1, 2, and 3 (Table|I]). Black 
lines: orthorhombic phases. Green and blue lines: monoclinic, 
sites 1 and 3 (see Fig. [2]). For each structure we calculated 
the occupations down the temperature at which the orbital 
polarizations are essentially complete {T ~ 200 K) and then 
extrapolated in a standard way 21] to T=0 K. 



(JT) elongation of a VO bond, that along y in cells 1 
and 4, and along x in cells 3 and 2, deforms the VOe 
octahedra. At about 140 in LaV03 and 116 K in YVO3, 
antiferroma^neizc C-type order develops (FM stacking of 
AFM a6-layers). At 77 K, YVO3 recovers the orthorhom- 
bic structure and the magnetic order changes from C- to 
G-type (3D- AFM), while the long VO bond becomes that 
along X in cells 1 and 3, and along y in 2 and 4. 

It has been suggested [l[ that these phase transitions 
are driven by the changes in a static orbital order (00) 
following the observed pattern of JT-distortions [1, Q ■ 
According to this JT-00 model, which assumes that the 
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FIG. 2: Primitive cell containing four ABO3 units. A-ions 
are orange, B-ions are at the centers of yellow O-octahedra. 
In terms of the primitive translation vectors, a, b, and c the 
global X, y, and z axes, directed approximately along the BO 
bonds, are: x = a/(2+2Q)+b/(2+2/3) y = -a/(2+2a)+b/(2+ 
2/?) and z = (c— c-b) / (2 + 27). Here, a, /3, and 7 are small and 
|x| = |y| = |zl =3.92A (3.82A) for LaVOa (YVO3). The B- 
containing fee-plane glide- mirrors (with translation (b — a) /2) 
unit 1 in 2 and 3 in 4, and exchanges the local, B-centered 
X and y coordinates. In the orthorhombic Pbnm structure 
(but not in the monoclinic P2i/a structure), the A-containing 
ab-plane mirrors unit 1 in 3 and 2 in 4. The empty crystal- 
field orbitals |3)^ of the monoclinic phase were put on sites 
1 = 1,2,3,4. Red (blue) indicate positive (negative) lobes. 
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Pbnm YVO3 (300 K) 
.56 -.21 .80 .96 
.83 .17 -.54 .53 
-.02 .96 .27 .51 
P2i/a YVO3 (100 K) 
.78 -.30 .55 .97 
.49 -.25 -.83 .58 
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Pbnm YVO3 (65 K) 
.64 -.26 .72 .98 
.72 .53 -.45 .60 
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TABLE I: LDA crystal-field (CF) levels wrt the t^g Fermi 
level, e^/meV (j = l,2,3), LDA CF orbitals at site i, in 
terms of the cubic orbitals, \xy)^ and \yz) in the global x, 
y, z axes defined in Fig. [5] Uj are LDA occupations. Orbitals 
at equivalent sites (see Fig. [2)|: (\f)A) is |ji)i (|j)3) with 
X ^ y; for the Pbnm structures is with z —z. 



crystal-field (CF) is due to the oxygen octahedra, the t2g 
orbital which is most antibonding with the O 2p orbitals, 
i.e. \sz) where s is the direction of the short in-plane VO 
bond, is empty; the other two t2g orbitals, due to Hund's 
rule coupling, are singly occupied. This 00 is C-type in 
the orthorhombic structure and G-type in the monoclinic 
structure. Later, this JT-00 model was challenged by a 
theory which assumes that the two highest orbitals, \xz) 
and \yz) , are basically degenerate so that orbital fluctu- 
ations play a key role [2, @]. Recently, ab initio LDA-I-C/ 
calculations gave support to the JT-00 model. 

In this Letter we show that in LaVOs quantum effects 
are strong down to 300 K, however they become negligible 
in the AF-C monoclinic phase. For YVO3 orbital fluctua- 
tions are suppressed already at 300 K, and the 77 K mag- 
netic transition is associated with a change of 00. We 
show that the CF splittings result not only from the JT-, 
but also from the GdFeOs-type distortions, and thus 00 
is intermediate between C- and G-type. The influence of 
the JT- and the GdFeOa-type distortion is, respectively, 
stronger and weaker than in the t\g titanates 13, 

The electronic structure is calculated with the 
LDA-I-DMFT method [igI fully accounting for the or- 
bital degrees of freedom [ll|]. First, we compute the LDA 
bands with the iV'^'-order mufhn-tin-orbital (NMTO) 
method [l7|; we obtain (for all structures) |-filled t2g 



bands separated by a ~0.5eV gap from the empty eg 
bands and by a ~2eV gap from the filled 2p bands. 
Next, we Lowdin downfold to Vi2g and remove the 
energy dependence of the downfolded orbitals by "iV- 
ization" 17|. These orbitals are strongly localized, hav- 
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ing Vt2g character only in their heads. Symmetric or- 
thonormalization finally yields localized [18| t2g Wannier 
functions and their corresponding Hamiltonian, H^^^. 
The many-body Hamiltonian is then a material-specific 
t2g Hubbard model, H = H^^^ + U ^ where for the on-site 
Coulomb repulsion, ?7, we use the conventional expres- 

U — 2 ^-^imcr m'cr' ^7na,7n' a' ^ima^im' a' t whcrC 

, and c|^^ creates an electron with spin a 
in a t2g Wannier orbital m at site i. The screened on-site 
Coulomb interaction is Uma,m'a' =USm,m'^<7,-a-' + {U' — 
JSa.<j'){^ — 5m,m'), whcrc J is the exchange term and 
U' — U — 2J the average Coulomb repulsion. We solve 
H in dynamical mean-field theory (DMFT) [2^, using a 
quantum Monte Carlo 21| impurity solver and working 
with the full self-energy matrix, "Smm' (w) Note that 
inversion is the only point symmetry of the V sites. The 
spectral matrix on the real w-axis is obtained by analytic 
continuation We use U =5 eV and J=0.68 eV, values 
close to theoretical [l^ and experimental [2^ estimates, 
which also give the correct mass renormalizations/Mott 
gaps for orthorhombic tl„ V/Ti oxides using the same 



computational scheme IG. Ill| 



Let us start by describing the LDA t2g bands in the 
orthorhombic 300 K phase. Remarkably, the CF orbitals 
\j)i (Table m j = 1,2,3), obtained by diagonalizing the 
on-site i block of iJ^^^, the hopping integrals t*'*/ (Ta- 
ble [III, the t2g band-shapes and band- width W (Fig. [3|) 
are rather similar to those of the t2g titanates. These 



similarities [24| are due to the similarity of the crystal 
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structures. Like in the titanates, the CFs are essentially 
determined by the GdFeOs-type distortion, mainly via 
the A-ligand field, specifically the AB and AOB cova- 
lency. However, in the vanadates the CF splittings are 
about half those of the respective titanates and the CF 
orbitals [25| and the hopping integrals are less deformed 
by cation covalency 2J]. This is due to the Ti — > V sub- 
stitution [2^ : since V is on the right of Ti in the periodic 
table, the V 3d level is closer to the O 2p and further from 
the A d level that the Ti 3d. Thus the sensitivity of the 
Bt2g Wannier functions to GdFeOa-type distortions de- 
creases, while the sensitivity to JT increases. 

Now, turning on the Coulomb repulsion transforms the 
metallic LDA density of states (DOS) into the spectral 
matrix of a Mott insulator (Fig. [3]). For LaVOs, the Mott 
gap is ~ 1 eV, in accord with optical conductivity data 
27|, and the Hubbard bands are centered around -1.5 eV 
and 2.5 eV, in very good agreement with photoemission 
and inverse photoemission For YVO3, the gap is 

slightly larger, ~1.2 eV, in accord with optical data [27|, 
and the Hubbard bands are centered around -1.5 eV and 



3 eV, in agreement with photoemission [28 1. 

The Mott gaps in the vanadates are somewhat larger 
than in the titanates, for which the measured gaps are 
-0.2 eV in LaTiOa and^ 1 eV in YTiOa jlj, in hne with 



LDA-I-DMFT results [ll|. This could appear surprising: 
orbital degeneracy increases the critical ratio for the Mott 
transition, Uc/W, by enhancing the effective band-width, 
and the enhancement is stronger the closer the system is 
to half- filling [2§|. So the gap should be smaller for a t^g 
than for a system, everything else remaining the same. 
However, the Hund's rule exchange energy, J, strongly 
suppresses this enhancement, as shown for half- filling in 
Ref . fsol. For n=^ and n=^, and using a 3- fold degenerate 
Hubbard model with a rectangular DOS, T=770K, and 
J/l¥~i (like in the vanadates where J/M^=0.68/1.9), 
we find that the metal to insulator transition occurs for 
U'/W ~ 1.5 when n=i, and for U' /W ~ 1.3 when n=i. 
So the Hund's-rule coupling dominates, and thus the van- 
dates can have larger gaps than the titanates. 

Like for the titanates [ll[, diagonalization of the 
DMFT occupation matrix yields three eigenvectors 
nearly identical to the LDA CF orbitals. For LaVOs 
at 770 K, the Coulomb repulsion only slightly increases 
the orbital polarization by changing the occupations as 
follows: 0.78^0.87, 0.63^0.65, and 0.59^0.48. Thus, 
surprisingly, orbital fluctuations are sizable and remain 
so down to room temperature: 713=0.26 at 290 K. Due to 
the stronger cation covalency in YVO3 , the Coulomb re- 
pulsion causes substantial orbital polarization already at 
770 K (see Fig.[l]). At 300 K, we find that only 4^4^ |0), 
paramagnetic with S=l, is occupied. Thus, YVO3 is or- 
bitally ordered well above the magnetic phase transition; 
since, at site 1, |3}i « \xz) (see Table|I| is empty and thus 
\xy) and \yz) are « singly occupied, the 00 happens to 
agree with the prediction of the JT-00 model [4. [gI. [i3| 
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FIG. 3: (color on-line) LDA+DMFT spectral matrix 

in the crystal-field basis. The off-diagonal terms are ~ 5 times 

smaller than in LDA. In black, the LDA DOS. 



(|3)i «-|3)3 ^\xz), 13)2 w-|3)4 « \yz)), even though 
the CF is caused mainly by the GdFe03-type distortion. 

What happens in the JT-distorted low temperature 
phases? The rms values of the hopping integrals hardly 
change, so W remains ~ 1.9 eV, but individual hopping 
integrals do change, even in the \xy) , \xz) , \yz) represen- 
tation. Most affected are the CF orbitals (see Table |l]). 

For LaV03, the CF splittings increase; this, in addition 
to the low temperature, lets the Coulomb repulsion sup- 
press quantum effects entirely. At sites 1 and 2, the oc- 
cupied state is in accord with the JT distortion (at site 1 
\xy) and \yz) are singly occupied and |3)j^ ~ \xz) is empty; 
at site 2, by symmetry, the empty state is [3)2 ~ \yz)', see 
Fig. [2] and Tab. [l| but this is not the case at sites 3 and 
4 (13)3^— \yz), 13)47^— \xz)). A static mean-field calcu- 
lation (pseudopotential-based LDA-|-t/) ^ yields empty 
states not far from ours [31], but, without analyzing the 
results, the 00 was ascribed to the JT distortions. In 
contrast, we find that the CF orbitals depend crucially 
also on the GdFe03-type distortions and that the 00 is 
not of G-type, but is intermediate between C- and G- 
type. 

For YVO3, the CF splittings are similar to those of the 
300 K phase, but quantum effects are negligible (Fig.[T]). 
On all sites in the monoclinic structure the empty orbital 
is almost the same as in the orthorhombic 300 K phase so 
00 does not follow the JT distortions {\3)^ « — [3)3 « 
\xz) , 13)9 ~ — 13)4 ~ \yz)), but is almost C-type. In the 
orthorhombic 77 K phase, the empty orbital at site 1, 
|3)j^ , only roughly equals \xz) . Our results are consistent 
with resonant x-ray scattering [5] and magnetization [l[ 
data. LDA-I-J7 yields results close to ours (sij. 

Finally, we find that the monoclinic structure favors 
C-type magnetic order over G-type by increasing some 
hopping integrals t^'*, (Table ITT|) to the empty orbital |3) 
along the c-direction. Assuming complete 00, conven- 
tional theory yields, for the superexchange couplings. 
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with j,j' CF orbitals and neighboring sites. We find 
that C-type order (./g^^^ < 0, J^^^ ■ 
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LaPbnm La P2i/a, site 1 La P2i/a, site 3 Y Pbnm Y P2i/a, site 1 Y P2i/a, site 3 Y Pbnm (65 K) 
j,j'\lmn 001 100 010 001 100 010 001 100 010 001 100 010 001 100 010 001 100 010 001 100 010 
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TABLE II: Hopping integrals t"y/meV from site i to a site i' — i + lx + my + nz, in the basis of crystal-field orbitals. 

Here i = l and (P2i/a only) i = 3. Notice that t':"+^ = and = For Pbnm structures only: = 

^ ' ■' ' 3,3 3 ,3 3,3 3 ,3 3,3 3 ,3 



vored over G-type, for which all couplings are positive, 
if J/U ^ 0.16. While the actual values of Jgg are sen- 
sitive to details 0, 0, this provides a microscopic 
explanation of C-type order in monoclinic LaVOa and 
YVO3, the change from C- to G-type across the struc- 
tural phase transition in YVO3, and thus could also ex- 
plain the magnetization-reversal phenomena [l[ . 

In conclusion, we find that the orthorhombic LaVOs is 
one of the few Mott insulators which exhibits large quan- 
tum effects at room temperature. This is not the case for 
YVO3 (and tig titanates [ll|). In the low temperature 
phases, orbital fluctuations are iiegligible for both vana- 
dates. This supports the view [ll.[3.ll4| that the magnetic 
structures of the vanadates can be explained by orbital- 
order. Recent LDA+U ^ and LDA-fPIRG calcu- 
lations agree with this, but previous literature ascribed 
00 mainly to JT-distortions. In contrast, we proved that 
both the JT and the GdFeOs-type distortions are crucial 
for the CF orbitals and their hopping integrals, and thus 
for the type of orbital and magnetic order. The effects of 
the GdFeO-type distortions are weaker and those of JT 
stronger than in titanates; their interplay is respon- 
sible for the rich phase diagram of the vanadates. 
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